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o : 

, We investigate some aspects of quintessence models with a non-minimally coupled scalar field 

■ and in particular we show that it can behave as a component of matter with —3 Pj p ^ 0. We 

' study the properties of gravitational waves in this class of models and discuss their energy spectrum 

and the cosmic microwave background anisotropies they induce. We also show that gravitational 
waves are damped by the anisotropic stress of the radiation and that their energy spectrum may 
help to distinguish between inverse power law potential and supergravity motivated potential. We 
finish by a discussion on the constraints arising from their density parameter f^Gw- 
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I. INTRODUCTION 



Recent astrophysical and cosmological observations such as the luminosity distance-redshift relation for the su- 
O ^ pernovae type la the recent observations of the cosmic microwave background temperature anisotropies [Q, 

gravitional lensing gand velocity fields |^ tend to indicate that a large fraction of the matter of the universe today 
is composed of matter with negative pressure (see e.g. for a comparison of the different observations). Recent 
analyses p,p| seem to indicate that the energy density p and the pressure P of this fluid satisfies 



- 1 < -P/p < -0.6, (1) 



k>( \ which is compatible with a cosmological constant A for which Pj p ~ —1 (see also ||T^ for arguments in favour 
. of Pj p < — 1). A typical value of J^a — 0.7 for its energy density in units of the critical density of the universe 
' corresponds to an energy scale of order 5 x 10~^^ GeV^ which is very far from what is expected from high energy 
physics; this is the well known cosmological constant problem To circumvent this problem different solutions 

have been proposed starting from the idea of a dynamical cosmological constant to lead to the class of models 
known as quintessence p^ , where a spatially homogeneous scalar field <f> is rolling down a potential decreasing when 
(\) tends to infinity. An example of such a potential which has been widely studied is the inverse power law potential. 
It can be obtained from some high energy physics models, e.g. where supersymmetry is broken through fermion 
condensates Q. Recently, it has been argued that supergravity has to be taken into account since today one 
expects the scalar field to be of order of the Planck mass Mpi and corrections to the potential appear at this energy. 
This leads to a better agreement with observations 

An important point about this family of models is the existence of scaling solutions |p7| , |l8| (refered to as tracking 
solutions) , i.e. such that (j) evolves as the scale factor of the universe at a given power. These solutions are attractors 
of the dynamical system describing the evolution of the scale factor and of the scalar field. This implies that the 
present time behaviour of the field is almost independent of its initial conditions This property allows to 

address |pT[| (i) the coincidence problem, i.e. the fact that (f> starts to dominate today and (ii) the fine tuning problem, 
i.e. the fact that one does not have to fine tune the initial condition of the field 4>. 

One of us extended these models to include a non-minimal coupling £^Rf{(f>) between the scalar field and the scalar 
curvature R ||2^]. Such a coupling term appears e.g. when quantising fields in curved spacetime p3|,p^ and in multi- 
dimensional theories [p5|-p7|. It was shown that when /(^) = ^^/2 tracking solutions still exist p2[~and this result 
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was generalised |28 to any coupling function / and potential V satisfying cx /"(^. However, such a coupling 

is constrained by the variations of the constants of nature |^9| which fix bounds on ^ |3C[|. A way to circumvent this 
problem is to consider quintessence models in the framework of scalar-tensor theories | |3l| , |32| where a double attractor 
mechanism can occur, i.e. of the scalar-tensor theory towards general relativity and of the scalar field (j) towards its 
tracking solution. 

Among all the possible observations of cosmology, gravitational waves give an insight on epochs where there was a 
variation of the background dynamics since every such variation affects the shape of the stochastic graviton background 
spectrum . We can then view our universe as containing a sea of stochastic gravitational waves from primordial 

origin, as predicted by most models of structure formations such as inflation ]33y34| ] (see also for a review) and 
topological defects scenarios [Q. Their spectrum extends typically from 10"^Hz (for wavelengths of order of the 
size of the Hubble radius today) to ~ 10^° Hz (the smallest mode that has been inflated out of the Hubble radius) 
and they could be detected or constrained by coming experiments such as LIGO VIRGO (at ~ 10^ Hz) 
and LISA (at ~ 10"'' Hz). Gravitational waves, which are perturbations in the metric of the universe have also 
an effect on the cosmic microwave (CMB) temperature anisotropy and polarisation |45| allowing to extract 

information on their amplitude from the measure of the CMB anisotropics. For instance, bounds on the energy density 
spectrum of these cosmological gravitational waves in units of the critical density, ilcw, have been obtained from the 
CMB [Eoirai 
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Gravitational waves are also a very good probe of the conditions in the early universe since they decou ple early in 
its history and can help e.g. testing the initial conditions of (j). An example was put forward by Giovannini ||46| , |47|| who 
showed that in a class of quintessential inflation models there was an era dominated by the scalar field (/> before 
the radiation dominated era which implies that a large part of the gravitational wave energy of order ficw — 10"^ 
(about eight orders of magnitude higher than for standard inflation) was in the GHz region. This may happen in 
any scenario where the inflation ends with a kinetic phase ||2^ , ^9| or when the dominant energy condition is violated 
pO|. On the other hand, the CMB temperature fluctuations give information on the history of the gravitational waves 
in between the last scattering surface and today through the integrated Sachs- Wolfe effect, whereas the polarisation 
of the CMB radiation gives mainly information on the gravitational waves at decoupling. These three observables 
(energy spectrum, CMB temperature and polarisation anisotropics) are thus complementary and we aim to present 
here a global study of the cosmological properties of the gravitational waves. 

The goals of this article are (i) to study in more details the cosmology with a non-minimal quintessence field and 
(ii) to study gravitational waves in this class of models. In we set up the general framework and describe the two 



potentials we shall consider. In §111 we introduce and define the observable quantities associated with the gravitational 
waves: their energy density spectrum and their imprint on the CMB radiation anisotropics and its polarisation. In 
IV , we point out the general mechanism of damping by the anisotropic stress of the radiation. In §^ we discuss the 



parameters of the problem and investigate the tuning of the potential parameters; we also describe the evolution of the 
bac kgro und spacetime and show that a non-minimally coupled quintessence field i s a ca ndidate for a (oj < — l)-matter. 



In §VH we describe the main properties of the gravitational waves. We finish in §VIII by presenting numerical results 
and we underline the complementarity of the different observational quantities. 

This work gives a detailed study of the observational effects of gravitational waves in the framework of quintessence, 
including some recent developments, and allowing for non-minimal coupling. This extends the work on quintessential 
inflation by including the effects on the CMB. It also extends the studies on ACDM |5|] to quintessence and is, 
as far as we know, a more complete study of the effect of gravitational waves on the CMB polarisation. We hope to 
show that a joint study of the gravitational wave detection experiments [^-|39| of the CMB experiments |^,^,^| and 
of the polarisation experiments can lead to a better determination of their properties. 

II. GENERAL FRAMEWORK 
A. Background spacetime 

We consider a universe described by a Friedmann-Lemaitre model with Euclidean spatial sections so that the metric 
takes the form 

ds^ = a^{r]) [-dr]"^ + S.jdx'dx^ = g^Ax^dx", (2) 
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where a is the scale factor and rj the conformal time. Greek indices run from to 3 and latin indices from 1 to 3. 

We assume that the matter content of the universe can be described by a mixure of matter and radiation (mainly 
baryons, CDM, photons and three families of massless, non-degenerate neutrinos) and a scalar field (p non-minimally 
coupled to gravity evolving in a potential ¥{((>) that will be described later. The action for this system is 



S 



R 



1 



ZK z 



(3) 



with K = SnG, G being the Newton constant, and where >Cmatter is the Lagrangian of the ordinary matter which is 
uncoupled to the scalar field and f{<j)) is an arbitrary function of the scalar field that will be specified later. The 
action (||) can be rewritten under the interesting form 



S = 



R 



with 



(4) 



(5) 



The stress-energy tensor of the scalar field is obtained by varying its Lagrangian [—£,Rf{4>) — ^df^(pd^<j) — V{<j))] to get 
Tl^f = ^f.4>^uc^ - ^ffp.VA</)VV - V{(j))g^j,, + 2C [5,..Va</)VV - V^0V.0 - 0V,,V,(/) + (j^Ucjjg^, + G^,/(0)] (6) 

where G^^, is the Einstein tensor of the metric g^^, V its covariant derivative and □ = V^V'^. 

The equations governing the evolution of the background spacetime are then obtained by varying (|^) with respect 
to gfiui 4> f-nd the ordinary matter fields to get respectively the Friedmann equations, the Klein-Gordon equation and 
the fluid conservation equation 
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(8) 
(9) 
(10) 



A dot denotes a derivative with respect to the conformal time and 7i = a/ a is the comoving Hubble parameter. The 
matter fluid energy density p and pressure P are assumed to satisfy the equation of state P — top. The factor oj varies 
from 1/3 deep in the radiation era to in the matter era. The scalar field energy density p^ and pressure P^ are 
obtained from its stress-energy tensor m) and are then explicitely given by 



P^ 



1^ 
2^ 

2a2 



V[<t>) 

v{4>) 



a? 

2e 



(27^ + 37^2)/(0)+7^/((/)) + /(0) 



We stress that the conservation equation derived from (^) reduces to the Klein-Gordon equation 
component, X say, we introduce the density parameter U,x defined as 



X 



KO^PX 

3^2 ■ 



(11) 

(12) 

For each matter 
(13) 



To completely specify the model, we have to fix the potential V^(0). Following our previous work and as 
discussed in the introduction we choose it to behave as 



Vi^) = AM - 



A 



a > 0. 



(14) 
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where A is an energy scale. As shown in [|2[ , such a potential leads to the existence of tracking solutions whatever the 
value of ^ and for which the scalar field behaves as a barotropic fluid of equation of state (as long as the background 
fluid dominates) 

P^ = uj^p^ with ^0 = -IH ^j^2 ' 

We also consider another class of potentials arising when one takes supergravity into account and given by 

= A-* exp(K?iV2) , a > 0. (16) 

The effect of the exponential term is important only at late time so that the scaling properties of the tracking solution 
are not affected during the matter and radiation era. However when the field starts to dominate it leads to its 
stabilisation which has an effect on the effective equation of state of the scalar fluid. 

B. Gravitational waves 

In this article, we want to focus on the properties of the gravitational waves which are tensorial perturbations. At 
linear order, the metric is expanded as 

9tiu = Qt^u + f,iu (17) 
where /^^ is a transverse traceless (TT) perturbation, i.e. satisfying 

/oo = /o. = o, f^.ur"^^, v^r"-o. (18) 

It is also useful to define the perturbation /i^j, by 

/^^ = a^h^^ (19) 

which, from (p^, satisfies 

/ioo = /»o. =0, /ifci<5'=' = 0, afc/i'='=0. (20) 

The equation of evolution of hi^i is obtained by considering the TT part of the perturbed Einstein equation (see 
e.g. |M) which leads to 



Ahki = 2KPa^TrM (21) 



where A = 9^9' is the Laplacian and where the anisotropic stress tensor of the matter, tt^;, is defined as the tensor 
component of the matter stress-energy tensor 

ST;^Pn^, nt^dXk^O. (22) 

The anisotropic stress of the matter fluid is dominated by the contribution of the neutrinos and of the photons and 



its form can be obtained by describing these relativistic fluids by a Boltzmann equation |^5| (sec § 111 B below). 



III. OBSERVATIONAL QUANTITIES 

The goal of this section is to deflne the observable quantities related to the gravitational waves. We start by 
reviewing the computation of the energy density pcw of a stochastic background of gravitational waves and finish 
by describing their effect on the CMB radiation, namely we present the computation of the coefficients Ci of the 
development of the angular correlation function of the CMB temperature anisotropy and polarisation. 



4 



A. Gravitational waves energy density 

The definition of tlie gravitational waves stress-energy tensor i^i^ can be found in e.g. |pl]] in the case of a Minkowski 
background spacetime, in e.g. ||5^ in the case of a Friedmann-Lemaitre spacetime and a general discussion can be 
found in e.g. p7[ . 

To define the gravitational waves stress-energy tensor, we have to expand the Einstein-Hilbert action (^) to second 
order in the perturbation /^j^, which implies to develop the curvature scalar R at second order in the perturbations 
(see [ pT| , ^6|j57| ) to get (up to divergence terms and forgetting the contribution arising from £mattcr) 

^^'^^ = - / -r\-yjo.0V''r^V^d'^- (23) 

J 4:Kcsm 

This expression is valid whatever the background metric as long as /^i^ is a transverse traceless perturbation. Note 
that contrarily to the "standard" situation, k now depends on (j) because of the non-minimal coupling with the scalar 
field. Using the fact that V^/q^ = o?dfj^hap we can rewrite the previous expression as 

5^^^S=- f -l—d^hkid^h'^'^d^x, (24) 



which assumes a Friedmann-Lemaitre background and the decomposition (pO[). Now, we decompose h^i on its two 
polarisations as 

hu^ E h(^Hfj,K)ei^\K), (25) 

A=+,x 

where e^.'J'^ (x) is the polarisation tensor defined as 

(x) ^ {elej - elef) S', + {def + ejel) (26) 

for a wave propagating along the direction 63 and where (ei, 62, 63) is a local orthonormal basis. Since this basis and 
the polarisation tensor satisfy 

^^^1 = 51 ei^\li,^=2S'y, (27) 
we can rewrite the action (|4|) of the graviton as 

- - E / Tr\nd,h('^d^hW^d^:,, (28) 
J 2kcs[(P\ 

which is the action for two massless scalar fields h\ evolving in the background spacetime, as first noticed by Gr- 
ishchuk p3| , ^ . By varying this action with respect to the background metric, we then deduce the stress-energy tensor 
of the gravitational waves 

= -^^AtT E {^^.h^^^^,h(^^ - g^M^^^d^h^^^) . (29) 
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If we decompose h'^^^ in Fourier modes as 

^/.(^)(k,r;)e^''-, (30) 

we can relate /i*^'*'-'(k, 77) to its initial value W-^^^iVi, rjiri) , i-e. its value deep in the radiation era (e.g. at the end of the 
inflationary phase) through the transfer function T{k,r]) by solving (21) to get 

/i(^'(k,r?) -T(fc,r;)/i(^)(k,??in). (31) 
Defining the initial power spectrum of the tensor modes as 
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{6 is the Dirac distribution), we can express the space average of io(x, rj) as 



2Koff [(/)]a2 ^ \ / Kcff[(/'] J 2Tra 

where we used an ergodic hypothesis to replace the space average by an ensemble average. Now, since (tQ(x, 77)) 
oscillates, we define the energy density of the gravitational waves by taking the average of (|33| ) over n periods. It 
follows that 

PGwiv) = -\n [ ^n(fc)T2(fc,,y)dln(fc), (34) 

where T(k, rf) is the root mean square of T(A;, rj) over n periods which is well defined as long as the amplitude of the 
wave varies slowly with respect to its period. 

The energy density pcw and energy density parameter ficw by frequency band are then obtained (after averaging 
on several periods of the wave) by 



dpGw(fc,?7) 1 



dln(fc) 27r2Keff[(/)] 



2 

2 



Ph{k)T'{k,7^), (35) 



dflGwik,!]) if K \ f k 



2 



din(fc) -en^wmjyn) ^''(^)^^(^'^)- (36) 

Let us stress some important points. Since we have to average on several periods, these expressions are valid 
only in a "shortwave limit" (see for discussion) where (i) the amplitude of the perturbation is small and (ii) the 
wavelength of the wave is small compared to the typical radius of the background spacetime. In our case this can be 
rephrased as k/Ti. > 1 which implies that the expressions ( p5[ - |3^ ) are valid only for modes which are "subhorizon" 
today, i.e. whose wavelength is smaller than the Hubble radius today. For such modes the ergodic hypothesis is well 
justified. In fact, because of the averaging procedure of the transfer function, wc have to restrict to modes such that 
k/Ho ^ 60 if we want to average on about ten periods. Again, we emphasize that there is an explicit dependence of 
f^GW E^nd pgw on the scalar field cj) because of the non-minimal coupling and our expressions reduce to the standard 
ones pl] , ^ , ^ , [l4|j4(jj57[| when ^ = 0. We have described the gravitational waves by two stochastic variables /i^'^^ 
which can be understood as being the classical limit of a complete quantum description of the gravitational waves 
(see e.g. ||,||,|6l for details). 

Before turning to the effects of the gravitational waves on the CMB, let us make a comment that will lead us 
to introduce some new notations. In the previous analysis we decomposed the metric perturbation hij on the basis 
Q^j{:s., k) = e^jexp(ik.x) of TT eigenfunctions of the Laplacian, i.e. such that 

AQ^j=~k'Q^^ with d'Q^^^S'^Q^j=0. (37) 

Such a decomposition is indeed not unique and we could have chosen any other such basis. In the CMB literature, 
one often prefers to use the basis Q^^^^ (x, k) defined by 

Qlf ) = -^1(61 ± ie2),(ei ± ie2),^'^-^, (38) 
the vectors ei and 62 being defined above (^. If we decompose hij on the latter basis as 

^31 



m=±2 ' 

then the two decompositions are related by 

3 



E / (|;^2i7(™)g(7)(x,k), (39) 



The two polarisations _ff ^^2) are then solution of (pl^) which reads 



(40) 
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(41) 



where the coefficient of the development of 5Tij as in (|39|) so that the transfer functions for and h''^^ are 

the same. If we now define the power spectrmxi of -ff (k, 77i„) as 



ff(™^)(k,r;in)i?^"^)*(k', 



'7in, 



{2TTyk-'PT{k)5{\i-\^)5m„ 



(42) 



one can easily check that 



Prik) = -Ph{k) 



(43) 



and that if the two polarisations + and x are independent then and ij' ^' are also independent. With these 

notations, the energy density spectra are given as 



dln(fc) 

dr^Gw 

dln(fc) 



1 



PT{k)f\k,l^), 



27r2 \k. 



= TT — hu PT{k)T\k,^) 



n 



(44) 
(45) 



Indeed, this does not change the result but we found interesting to make the link between the notations used in 
the gravitational waves literature ^1 35 4^,^,^,^ and in the CMB literature |45|j5^ ] , specially because we want to 
present both in a unified framework and language. From now on, we use the second decomposition and its interest 
will be enhghtened by the study of the CMB anisotropies. 



B. CMB temperature and polarisation anisotropies 

Gravitational waves, being metric perturbations, have an effect on the temperature and polarisation of the CMB 
photons. Any metric perturbation induces a fiuctuation on the CMB temperature Q through the Sachs- Wolfe effect |5|] 
and any anisotropic distribution of photons scattered by electrons will become polarised and vice- versa. Since Thomson 
scattering generates linear polarisation, we only need to consider the Stokes parameters Q and U and more conveniently 
their two combinations Q ±iU which are invariant under rotation. 

Following Hu and White js^ , we decompose the tensorial part of the temperature anisotropies according to 

e(,7,x,n)= / E ef)(fc,,7)Gr(k,x,n), (46) 

{Q±^UM^,h)= I E (i?^^±*i3l"^)(fc,^)±2Gr(k,x,n). (47) 

The coefficients E^™''^ and B'f^^ transform as Eg {—YEg and Bg ~{~YBg under parity and are called the "elec- 
tric" and "magnetic" part of the polarisation. The functions G™, ±2G"' form three independent sets of orthonormal 
functions and depend both on the position x and on the direction of propagation of the photons n and are defined as 



Gr(k, X, n) ^ i^Ytin) exp(zk.x), (48) 



±2Gr(k,x,n) ^ (-)y-i^ ±2y,"(r^)exp(zk.x), (49) 

where the functions Y^{h) are the standard spherical harmonics and the functions ±2yi"{fi') are the spin-weighted 
spherical harmonics p0|-|62[. Note that the decomposition on the basis Q^f^^ in the previous section is enlightened 



by the fact that Q^J'^n'ni = G^"'^ i|,|5|. 

The angular correlation function of these temperature/polarisation anisotropies are observed on a 2-sphere around 
us and can be decomposed in Legendre polynomials Pi as 
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5C/. 
T 



(71)^(72)) - ^^(2£+ l)CrPK7i-72), 



(50) 



where (7, V stand for G, i5, or B. Now the brackets stand for an average on the sky, i.e. on all pairs (71,72) such 
that 71.72 = cos 012- Using the orthonormality properties of the eigenfunctions G, equations (|4^^) can be inverted 
to extract the angular power spectra C^^ of the temperature and polarisation anisotropies as 



•' m=±2 



(51) 



The equations of evolution of the temperature and polarisation multipoles i?^™-* and i?^™^ can be obtained by 

decomposing the Boltzmann equation satisfied by the photon (or neutrino) distribution function on the eigenfunctions 
G to get 



e 



k 



QKg (m) _ Ql^e+l ^(m) 



M""^ = k 



b'T^ = k 



2^-1 



Tef^ + '5.25("), 



2™ Irn) _ _2K^ (,„) 

2^ + 3 



2m („j) _ 2^^i (m) 
2^ + 3^+1 



where 



l{l + 1) 



5'('") = f p(™) _ ^(™) 
10 V 2 



(m) 



(m) 
'2 



1 - 



(52) 
(53) 
(54) 

(55) 
(56) 

(57) 



The differential optical depth r vanishes for neutrinos (except in the very early universe, but the observable wavelengths 
in the CMB are not affected by this) and is proportional to the free electron density in the case of photons. It has to 
be calculated by solving the relevant kinetic recombination equations for hydrogen and helium |]6^|65|| . 

The quantity p(™) represents the coupling between temperature and polarisation. Due to our choice of decomposi- 
tion, only the electric part of the polarisation is affected by temperature anisotropies. However, electric and magnetic 
part of polarisation couple themselves as photons propagate. The Clebsch-Gordan coefficients s>^T arise from product 
properties of spherical harmonics. They are obtained in the same way as for the scalar modes |65| . 



1)0 



ef'(r?o,fc) 

2^+1 

igf^(^o,fc) 
2^+1 

= /" d, re-F(™)/3('")(fc(,o - v))- 



dr;e--5(")j^'")(fc(ryo-r;)), 
^~V6 d?7 fe-"p("')e^'"^(fc(?7o - v)), 



(58) 
(59) 
(60) 



The functions j^'"'', e^™-* 



and Z?^"^ are defined in terms of the spherical Bessel functions, j{{x) as 



l3i£ + 2y.,n{x) 



8(^-2)! a;2 



(61) 
(62) 
(63) 



where a prime denotes a derivative with respect to the argument x. 
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Beside the small contribution due to the polarisation in Eq. (p8|), the temperature fluctuation in the direction 7 
reduce to the well known result by Sachs and Wolfe 

5T I r° ■ 

-7f,{l)^-7; /iy7Vdry, (64) 

^ ^ "''JLSS 

where the subscript LSS stands for last scattering surface. The "visibility function" fe~'^ appearing in equations 
( p9[ - |60| ) takes a non zero value only at the time of decoupling so that, contrarily to temperature anisotropics which are 
constantly generated by gravitational interactions with photons, polarisation is generated only at the last scattering 
surface. 



IV. DAMPING OF THE GRAVITATIONAL WAVES 



We first study the effect of the damping of the gravitational waves due to the anisotropic stress of the photons. 
Such a damping of the amplitude of gravitational waves in various viscous cosmic media has been already discussed 
in 1?^]; we give a description of this damping in the formalism we use here in order to quantify precisely its effect on 
CMB anisotropics. On subhorizon scales larger than the diffusion length Ad = of the photons, i.e. such that 

Hofi' EE T.^^ < fc < f, (65) 

the set of equations ( p2[ - p7| ) implies that 

er^=-^| and i^(±^)=-^ep). (66) 

Since tt'^^^-' is proportional to 02^^\ we can insert the former expressions in the gravitational waves evolution equa- 
tion (|4^) to get the back reaction of the anisotropic stress 

32 II 

H + 2ncsH + k^H= nPa^ — . (67) 

15 f 

Setting 

H = Hh, (68) 
with H solution of the homogeneous equation of evolution (i.e. with tt^™' =0), we get 

H'h+ \2H + 2rLcsH + —kPo^— ] h + —uPa^ — h = 0. (69) 
\ 15 f / 15 f 

Now, (i) since H ~ kH, we deduce that H HcsH, (ii) since nPa^ ~ Ti-cB, it follows that (32/15)KPa^ff/'r ~ 
Ti^ff/r <C kH ~ H and (iii) since h ~ h/TD, Hh ^ Hh and in conclusion in the limit (|6^) the equation of evolution 
of the gravitational waves in presence of the anisotropic stress (^9[) reduces to 

16 

h = KPa'^ih. (70) 

15 ^ ' 

We deduce that a mode k is damped from the time it enters the Hubble radius, i.e. ij ~ k~^ since it happens during 
the radiation era, to roughly the time when the anisotropic stress becomes negligible, i.e. approximately at the time 
of last scattering ?7lss- It follows that 

/ 1 a /■'7LSS ^p„2 \ 
hik, 77LSS) =^ exp ~- / —^dr, h (fc, 1/fc) . (71) 



15 



i/fc 



This damping of the gravitational waves by the anisotropic stress of the photon fluid is analogous to the damping of 
the scalar modes (density fluctuations) known as the Silk damping ]7^ a description of which, in the formalism used 
here, can be found in [p8|. Note however that the origin of the damping is different. 
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This effect is small but, apart from [Q, was not much emphasized in the hterature before. Assuming that the 
universe is completely ionized until the last scattering surface, the integral of Eq. (|7^) is of order |Q 



Ul - ^ ) 7^^«0(^LSS - l/k) 10-%1 1/fc^LSSj 



(72) 



where uip is the proton mass, UTh is Thomson scattering cross section. The real damping factor is greater than the 



estimate ( [72[ ) because the universe becomes neutral at the last scattering surface (so that the term f is smaller). 
In Fig. |l| [left], we plot this damping factor for the modes that entered into the Hubble radius long before the last 
scattering surface (i.e. such that k ^ '7ls^s)- ^ consequence, the comparison between the damped case to the 
undamped case, shozn on Fig. |^ [right] does not show significant differences. The amplitude of the high-£ tail of the 
CMB anisotropy spectrum is lowered by roughly 10% when one includes this effect. The same occurs of course for 
polarisation. We emphasize that this result does not depend on any particular model, and is not included in the most 
recent (3.2) version of CMBFAST. 
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FIG. 1. Damping of the gravitational waves due to its coupling to the photons anisotropic stress. The damping (left figure) is 
generated only when the universe becomes neutral, i.e. soon before the last scattering surface. As a consequence, all the mode 
which have already entered into the Hubble radius at recombination are equally damped, regardless of their wavelength. The 
influence of this damping on the CMB anisotropics is shown on the right flgure. Since all the modes are equally damped, this 
translates into a constant ratio of amplitude between the damped and the undamped cases. Note however that our derivation 
is valid only when Eq. ( |6^ ) applies, which is not true at the end of decouling, when r becomes small. This is why the actual 
damping (10%) is smaller than what is expected from the left plot. 



V. SPECIFICATION OF THE MODEL 



A. Model parameters 



At this stage, the model we are discussing depends on five parameters: 

1. /((/>) which is an arbitrary function of the scalar field (f), 

2. ^ the coupling of the scalar field with the background spacetime. 



3. a the slope of the potential (14) or ([16D, 

4. the energy density in the scalar field today, 

5. Prik) the spectrum of the gravitational waves. 

Indeed there exist some constraints on these functions and parameters and we make the following assumptions and 
choices: 

1. We assume that /(^) = ^^/2; this is the only choice for which the coupling constant ^ is dimensionless. Moreover 
such a choice can be seen as the lowest term in an expansion of / in powers of (p. As shown in p2| there exists 
tracking solutions for the field evolving in the potential dlj) with such a coupling. 
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If the scalar field (/) is coupled to the spacetime metric, this coupling must be weak enough so that it does not 
generate a significant time variation of the constants of nature |29[| . Taking into account the bound on the 
variation \G/G\ of the Newton constant and on the variation \a/a\ of the fine structure constant it 
was shown pO] that 

- 10^^ < C < 10"^ - 10^^ (73) 
This bound is however sensitive to the shape of the potential. On the other hand the experimental constraints 



(from the Shapiro effect and the light deflection in the Solar system) on the post-Newtonian parameters [g? 68 
imply ^ 

ICI < '4^. (74) 

However, in this class of models one does not try to have a theory converging towards general relativity at 
late time and the coupling ^ is constant which is the main reasons of the severe bounds on its value. This 
can be improved by generalising this kind of models by considering them in the framework of scalar-tensor 
theories 

3. In most models a is not constrained theoretically. If the matter content of the universe today is dominated by 
the matter-radiation fluid then the fact that the observations |^ favour —1 < ui^ < —0.6 gives a bound on a, 
which is indeeed not the case anymore if the scalar field starts to dominate. In Fig. ^ [left], we compare this 
analytic estimate and the numerical determination of the energy scale A as a function of the slope a. We see 
that if a > 4 then A is at least larger than 1 TeV (when fi^ = 0.7). 

4. The density parameter fi^ is not severely constrained theoretically, but observations seem to indicate fi^ ~ 0.7. 
One has to check that if the scalar field was dominating the matter content of the universe at some early stage 
then it has to be subdominant at the time of nucleosynthesis (see e.g. [^). The choice of fixes the value 
of the energy scale A in ( [l^ ) or (p^; this is the coincidence problem. On Fig. ^ [right], we depict the variation 
of the energy scale A with 57^ and a. It is not very sensitive to 17° as long as 0.1 < fi^ < 0.9. In fact, 
when the quintessence field starts to dominate the matter content and if we have reached the attractor then 
d?V/d(jP' (X (see [Q), and ~ V/M'^y so that we can estimate that the variation of A with a follows 

A=(pcritM^i)^- (75) 

We conclude that 

M 1 6^11 



A 4 + a 17^ 

and thus a prescision of 10% on il^ requires to tune A at a 1% level if e.g. a — 6, which is a less drastic tuning 
than the usual cosmological constant fine tuning problem. 




FIG. 2. Variation of the energy scale A of the potential ( |l4[ ) with the energy density of the scalar field fi^ and the slope of the 
potential a. We first show [left] the variation of A with a when £7^ — 0.7 and the comparison with the analytic estimate ( fz^ ) 
and [right] the variation of A with for a — 2, Q, 12. 
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5. Pt has to be determined by a specific model, such as e.g. inflation, and we parametrise it as 

Prik) = Ark"''^ (76) 

where At is a constant dans nx is the tensor mode spectral index. At is obtained by normalising the CMB 
temperature anisotropics to the COBE data at ^ = 10 for which 



Toy^^^CT^SO^K. (77) 

Since some measurements tend to show that there is a peak at the degree scale , we conclude that a significant 
part of the anisotropics may be generated by the scalar modes. In the "standard" slow-roll inflation picture, this 
is compatible with an almost scale-invariant spectrum with a low tensor contribution, in which case the COBE 
results would put only an upper limit on the amplitude of the gravitational waves spectrum. Nevertheless, we 
point out that it is also possible that most of the large scale anisotropics can be generated by gravitational waves. 
This assumes a strong deviation from scale invariance {ns = 1.69 and ut = 0.0), but is in good agreement with 
observational data IITCl . 



B. Initial conditions and behaviour of the background spacetime 

Concerning the initial conditions for the scalar field 0, we will consider the two extreme cases: 

• (ICl) where we assume that the scalar field is at equipartition with the matter (i.e. mainly with the radiation) 
deep in the radiation era, 

• (IC2) where it dominates the matter content of the universe at a very early stage. 
Situation (ICl) implies that at the end of reheating, 

P4> < (78) 

where the factor 10^* is roughly the inverse of the number of degrees of freedom at that time. Since the quintessence 
field is already subdominant at this epoch, one does not need to care about its effect on nucleosynthesis since it 
remains subdominant until recently. In the second situation (IC2), the field starts by dominating and inflation ends 
by a kinetic phase rapidly than and will thus become subdominant. One has to check that this happens before 
nucleosynthesis pqMS] . A realisation of such initial conditions can be obtained in quintessential inflation Ba] . 




1/(1 +z) 1/(1+z) 
FIG. 3. The behaviour of the energy density of the matter (long-dashed line), radiation (short-dashed line) and scalar 
field (solid and dotted lines), as a function of the redshift for the two class of initial conditions: (ICl), when the field is at 
equipartition with the radiation [left], and (IC2), when the field initially dominates the matter content of the universe [right]. 
The solid line represents the case when the field evolves in an inverse power-law potential, and the dotted line represents the 
case when the field evolves in the SUGRA potential. Note that when the field dominates at early times, the SUGRA potential 
stabilizes the field, which reaches the tracking solution earlier. 
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FIG. 4. The behaviour of the equation of state parameter as a function of the redshift for the two class of initial conditions: 
(ICl), when the field is at equipartition with the radiation [left], and (IC2), when the field initially dominates the matter 
content of the universe [right]. The solid line represents the case when the field evolves in an inverse power-law potential, 
and the dotted line represents the case when the field evolves in the SUGRA potential. In the case of (IC2), the field reaches 
the tracking solution only when SUGRA corrections to the potential are considered. Note also the spikes in the SUGRA case 
[right], which illustrate the fact that the field bounces around the Planck scale (see Fig. ^ below). 



10-= 



10" 



10"' 



S 10 



10-' 



-e- 



10"' 



10" 



10" 



10" 



Quintessence, 
SUGRA Quintessence, 



(101) 
(101) 



10^ 



10" 



10"' 



S 10 



10"' 



-e- 



10"' 



10" 



10" 



10" 



Quintessence, (102) 
SUGRA Quintessence, (102) 



10 



■30 



10 



25 



10 



■20 



.|Q-15 .|Q-IU ^Q-J ^QV ^Q-OV ^Q-^-J ^Q-^U ^ Q- 

1/(1+Z) 1/(1+Z) 

FIG. 5. The behaviour of the quintessence field as a function of the redshift for the two class of initial conditions: (ICl), 
when the field is at equipartition with the radiation [left], and (IC2), when the field initially dominates the matter content of 
the universe [right]. The solid line represents the case when the field evolves in an inverse power-law potential, and the dotted 
line represents the case when the field evolves in the SUGRA potential. In the case of (ICl), the field always reaches the 
tracking solution before today, whereas for (IC2), the field reaches the tracking solution only when it evolves in the SUGRA 
potential. 

In Fig. ^, we depict the evolution of the energy density of the quintessence field, matter and radiation for the 
initial conditions (ICl) [left] and (IC2) [right]. We see that for a very large range of initial conditions (roughly for 
lO-^^GeV* < p0 < IQi^GeV at a redshift of z ~ 10^°) we end up with a quintessence field which starts to dominate 
today. This explains briefly how the fine tuning problem is solved |2|]. We can also check that with these values 
the scalar field does not dominate the matter content of the universe at nucleosynthesis, i.e. at a redshift of order 
z~ IQio. 

An interesting point concerns the evolution of the scalar field equation of state in the case (IC2) when ^ ~ 0- The 
field rolls down very fastly so that we are first in a regime where 



10 



■30 



25 



■20 



,-15 



10 



10 



10"' 



10" 



p 



102 

2^ 



(79) 



13 



from which we conclude that its equation of state is Ci;0 ~ 1 (see Fig. ||). But, due to the exponential behaviour 
contribution of the potential, the field is stopped when (j) > Mpi and then rolls back to smaller values (see Fig. ^) so 
that the field undergoes a series of damped oscillations (because of the friction term coming from the expansion in the 
Klein-Gordon equation) . This implies that there exist times such that ~ and thus small period around them where 
the equation of state varies rapidly to w,/, ~ — 1 (see Fig. This sudden change in the equation of state of cf) happen 
while it is dominating the matter content of the universe (see Fig. ^) so that it implies variations in the evolution of 
the scale factor of the universe which, in principle, should let a signature in the gravitational waves energy spectrum. 
Indeed, this does not happen in standard quintessence and is a specific feature of the SUGRA-quintessence. 

When ^ ^ 0, there are no significant modifications to the background dynamics as long as the field has not reached 
the Planck mass [because 2^k/(0) is small compared to unity, see Eq. ^. Then, the main difference appears at late 
times when the field starts to dominate and comes from the fact that the bound — 1 < < 1 no longer applies, and 
one can get lower values of to^. Equivalently, the equation of state parameter to = (Pfluid + ^V)/(pfluid + P0) for the 
whole background fluids can reach values smaller than —1 (see Fig. where we plot the variation of oj^ as a function 
of redshift). As pointed out by Caldwell jl^, such a matter fits the current observational data. Different candidates 
such as a decaying dark matter component [ [zil and a kinetic quintessence field fz^ were proposed. Here, we show 
that any non- minimally coupled scalar field may be a good candidate for a component of matter with to < —1. The 
constraints (^3|) on ^ implies that for our class of models 

-3 < < 

if the scalar field dominates. We emphasize that uj^ is not uniquely defined according to the way one splits T^'^'^ in (|^). 
In Fig. we used the Friedmann equations (0j^) to extract oj from 

^ - 1 = --(l+Lo)n 

and then uj^ from 

i 

where i runs on all the matter species. This corresponds to the value of lo as it may be reconstructed from observational 
data such as e.g. the supernovae type la. 




FIG. 6. Evolution of the equation of state parameter u as a function of the redshift for different values of the coupling ^. As 
soon as the coup ling is not minimal, uj can reach values smaller than —1. The parameters of the model considered here are: 
a — &, potential (hq) including SUGRA corrections, — 0.7 and (ICl) initial conditions. 



VI. QUALITATIVE DISCUSSION 
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A. Gravitational waves spectrum 



Equation ( pl| ) describes the evolution of a damped oscillator. Injecting the ansatz 

jj(rn) ^ ^(m) exp(i/!c,y) (80) 



in (41) and performing a WKB approximation leads to the equation 

+ = (81) 

for the evolution of the amplitude where Hcs = a/a, with 

5=aVl-2</(0). (82) 

This WKB approximation holds only for "sub- horizon" modes. Before a mode has a wavelength smaller than the 
Hubble radius, its amplitude evolves according to 

+ 2Weffi^"^ = 0, (83) 

the solutions of which are a constant mode and a decaying mode. Neglecting the decaying mode, we see that the 
wave is "frozen" as long as its wavelength is larger than the Hubble radius, and that it undergoes damped oscillations 
once its wavelength is shorter than the Hubble radius. The damping of a mode of wavenumber k between the time it 
enters the Hubble radius and today is then proportional to 

^, (84) 
ao 

where is the scale factor evaluated at the time r]k when the mode k enters the Hubble radius (i.e. when Ti. — k) and 
ao is scale factor today. Injecting this behaviour in (|3^), we obtain that the energy density spectrum of gravitational 
waves scales as 

f^^<.eriPrik). (85) 

First let us assume that ^ = 0. For wavelengths corresponding to modes that have entered the Hubble radius in 
the matter dominated era (for which a oa i]^ and thus ry^ ~ k~^), one can easily sort out that 

ttk ^ k~^ (86) 

and the gravitational waves spectrum behaves as 

Equivalently, for wavelengths corresponding to modes entering the Hubble radius in the radiation dominated era (for 
which a (X rf) one can show that the gravitational waves energy spectrum behaves as 



dr^GW 
dln(fc) 



oc k^Prik). (88) 



To finish, if it happens that there exist wavelengths corresponding to modes that have entered the Hubble radius 
while the scalar field was dominating (for which a oc y/rj since oc 1/a^) one obtains that 



dr^GW 
dln(fc) 



oc 



k^Prik). (89) 



In conclusion, we have found three behaviours for the gravitational waves spectrum according to the wavelength. In 
Fig. 1^, we give an example of such a spectrum in a case where one has a scalar field dominating at early stage [initial 
condition (IC2)]. These results hold also when ^ but the slopes of the spectrum are slightly changed since the 
time behaviour of a has to be replaced by the one of a. 
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B. CMB anisotropies 



For scales smaller than the Hubble radius at decoupling, one can follow the same lines to predict the tensor part 
of the CMB temperature anisotropies. The main difference is that the expression for £{i + l)Ci does not involve any 
factor fc^ as in Eqns. (^), the reason being that Eq. can be integrated by parts to drop the time derivative of 

which shows that anisotropies are mostly generated on the last scattering surface with an amplitude of 
Therefore, the spectrum behaves as 

£(^ + l)Q(x^^-^^•^-^r^-^ (90) 

for modes which have entered the Hubble radius is the matter dominated, radiation dominated and kinetic scalar field 
dominated eras respectively. With standard cosmological parameters, the radiation to matter transition occurs soon 
before the decoupling, and the scalar field dominates only at very early times. As a consequence, one sees almost only 
the regime £{i + l)Cg oc For modes which enter into the Hubble radius after the last scattering surface, one 

can show that the produced spectrum scales as 

+ i)Q ocr^. (91) 

Note that this expression is indeed an approximation and that it is not easy to calculate an accurate analytical 
solution l?^. These results are illustrated on Fig. 0. As already stressed, the result of Eq. ^ applies at large angular 
scales which have not entered into the Hubble radius at recombination. For standard cosmologies, this occurs for 
multipoles smaller than t ~ 100 (in addition, there are also some corrections to this rough estimate which occur at the 
very smallest multipoles and slightly boost the spectrum, as can also be seen on Fig. Then, at higher multipoles the 
result of Eq. |9^ is valid. The matter dominated regime before recombination is rather short, and occurs only between 
I ~ 100 and i ~ 200 (less than one oscillation in the spectrum). For I > 200, one sees the regime i{i + \)Ci oc 
(see also Fig. 1 of Q). 



C. Results of the ACDM model 

Before turning to a more general numerical study of the class of models we consider in this article, we recall in Fig. |^ 
the general results for the temperature and polarisation angular power spectra and the gravitational waves density 
spectrum for a ACDM model. This spectrum has two branches: a soft branch at lower frequencies (corresponding to 
the matter dominated era) and a hard branch at higher frequencies (corresponding to modes that entered the horizon 
in the radiation era) . Following , we set the cut-off on this spectrum to the last mode that has been inflated out 
of the Hubble radius. 




10 100 1000 10"''^ 10"''° 10"^ 10° 10^ 10^° 
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FIG. 7. The temperature et polarisation of the CMB induced by gravitational waves [left] and their energy density spectrum 
[right] in a standard ACDM model with ut ~ 0. 



VII. NUMERICAL RESULTS 
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A. Field (j) initially at equipartition 



Since the scalar field only starts to dominate at very recent time, we expect no efi^ect on the gravitational waves 
energy spectrum (since at earlier time the universe is always radiation dominated). However, the change in today's 
universe equation of state yield a specific signature in the angular diameter-distance relation. Hence, one expects to 
see the quintessence field behaviour in the positions of the peaks in the CMB anisotropy spectra. 

The temperature anisotropics plots of Fig. |^ are therefore identical at high multipoles except for their overall 
position which are different. At low redshift, the scalar field dominates and the dynamics of the expansion depends 
explicitely of the value of the coupling ^, which cause some slight differences in the CMB anisotropics at the very first 
multipoles < 5). We have also seen that the polarisation is generated by gravity and therefore different gravitational 
constants lead to different normalisation between the polarisation and the temperature spectra. Since we normalise 
the "bare" Einstein constant k so that the effective Einstein constant corresponds to what we measure (in e.g. a 
Cavendish experiment), models with a different ^ have different n. At decoupling, the scalar field does not dominate 
and therefore i^^^ = k^^^. This induces different amplitudes for the polarisation anisotropy spectra. For the lowest 
values of ^ there is a factor 2 — 4 in amplitude as compared with the ^ = case, which roughly corresponds to the 
square of the variation of Kog (and, hence G) between the last scattering surface and now. Note that the effect of f 
depends on its sign. This is the reason why the constraint derived by Chiba are stronger for negative values of ^. 
The same can be seen in Fig. ^. 

We conclude that the temperature anisotropics and polarisation give mainly information on the spectral index tit, 
the energy density of the scalar field today and its coupling ^. 
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FIG. 8. Influence of the coupling ^ on the CMB temperature and polarisation anisotropies. 
angular diameter-distance relation and therefore affects the overall position of the spectrum. 



The value of ^ influences the 
The parameters of the model 



considered here are the same as in Fig. ^ a = 6, potential (^6|) including SUGRA corrections, 57^ — 0.7 and ut ~ 0. 



B. Field (j) dominates at early stage 

We now turn to the more unusual case where the scalar field dominates at the end of inflation and where the universe 
undergoes a kinetic phase before the radiation era (2^,^^ as in e.g. quintessential inflation The wavelengths 

corresponding to the observable CMB multipoles {t < 2000) are much larger than the Hubble radius at nucleosynthesis, 
epoch at which we have to be radiation dominated. As a consequence, we expect no signature from this early phase 
on the CMB anisotropies and polarisation. 

As first pointed out in |4^,^^, if the scalar field dominates at early stage, there is an excess of gravitational waves 
at high frequency [see equation (|89|)]. On Fig. ^ we present such a spectrum and we will discuss the implication of 
this excess later. 

An interesting effect concerns the difference between the spectra obtained from an inverse power law potential and 
a SUGRA-like potential. As shown on Fig. ^ [right], the amplitude of the spectrum at high frequency in roughly 
30% higher for inverse power law potentials. The relative decrease in amplitude at these frequencies for SUGRA-like 
potentials depends on the dynamics of the scalar field in the bounce (see Figs. ^ and H) during which the equation of 
state varies from -fl to —1 and to -f 1 again. Thus, during this time, the modes that had just entered into the Hubble 
radius (and thus which had just started to undergo damped oscillations) went out of it (during the uj < Q epoch) and 
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their amplitude was frozen before re-entering the Hubble radius again. Hence, the modes of larger wavelengths are 
less damped which explains this decrease in amplitude. Now, if the slope of the potential is less steep, the bounce 
lasts longer (note that we always reach = — 1 at the point where (p = 0) and thus the damping is stronger. This 
signature, even if not detectable by coming experiment is nevertheless a clear feature of supergravity. 

To finish, let us discuss the total energy density of gravitational waves f^Q^ today. As pointed out in , it has 
also to be negligible at nucleosynthesis; this constraint is more drastic than the only requirement that f2^ be negligible 
at that time. Let us emphasize that the constraint on fl^ cannot be avoided (since it involves background dynamics) 
whereas the one on f^Q^^ depends on At and ut and thus leads to a combined constraint on the initial conditions 
of the scalar field and on the initial power spectrum of the gravitational waves. Besides At and ut, ^q^v mainly 
depends on the initial values of and pi-ad which can be parametrised by the reheating temperature Tji (related 
roughly to Prad at that time) and the redshift of equality between the kinetic scalar field era and the radiation era 
(related roughly to p^/prad at the end of reheating). f^Qw '^^n be estimated by the surface of the spectrum below the 
part with a positive slope (i.e. the high frequency part; see Fig. ^) and thus of order 



(-,0 _ , dilcw 



(92) 



where and are respectively the modes entering the Hubble radius at the reheating and at z*. Thus the "bump" 
at short wavelength cannot be too high. Moreover, the energy density at the end of reheating cannot be higher that 
Planck scale, so that it fixes a limit on the shortest mode in which gravitational waves are produced. On Fig. |lO|, we 
first plot [left] the variation of the gravitational wave spectrum with the parameters (Tr, z») and we then give [right] 
the "safe" zone of parameters for nucleosynthesis [for ut = 0] and defined [Q,Q by ^Q^J^ ^ 10~^. Let us briefly 
explain how this bounds are obtained. 

1. We first rephrase in terms of Tr the fact that the field is dominating at the end of the inflation phase, i.e. 

z* < zr z* < aiTfl, (93) 

where ai is some numerical coefficient. This corresponds to the solid line on Fig. 

2. We then impose that the scalar field is subdominant at nucleosynthesis, i.e. that 

z* > 10^°. (94) 

This corresponds to the horizontal dash-dot line on Fig. ^ 

3. At the end of the inflation phase, we want the energy density to be smaller that the Planck energy density. If 
the scalar field is dominating if gives 

(1 + 1^]' (1 + ^eq) < PPi ^z*> a2Tl (95) 



where a2 is another numerical coefficient. This corresponds to the dot line on Fig. 10. Note that since we are 
in a field dominated era oc a? and thus on this "Planck limit" we have fc/j oc z'^Jz^, and A:* cx z* (and thus 
kR cx 1/zi? and fc* cx Zr) from wich we conclude that the maximum of the power spectrum is roughly located 
on a curve {1/zr, 1/z^j [see Fig. p^ . 

4. To finish, we want that the gravitational waves energy density does not alter nucleosynthesis, i.e. that 

f^GW ^ 10"' ^Tr< asz, (96) 

where is a third numerical coefficient. This corresponds to the dot-dash line on Fig. p^ . 

For all the points {Tr, z*) above the dotted and dot-dashed lines of Fig. [l^ [right] , there is no excess of gravitational 
waves. The solid line separates the two sets of initial conditions zwe have considered. We must emphasize that this 
result was obtained for ut = and that the spectrum can be tilted, which modifies the bounds on the parameter 
set (Tfj, z*) (more precisely, taking smaller At or negative ut lowers the diagonal dot-dashed line). Such constraints 
may be important for instance while considering models where a scalar field dominates at baryogenesis fl% . 

In the case of a "blue" initial power spectrum (i.e. with ut > in our notations, or /? > —2 in the notations 
of [|7^), as advocated for example in [Q, the flat branch of Fig. ^ (corresponding to the "semi-hard" branch of [|7|) 
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is already tilted, giving as stronger constraint on our model. For instance, if ut = 0.2, the amplitude at u = 10 GHz 
is boosted by a factor ~ 3 x 10^. As a consequence, the quantity of gravitational waves at high frequency cannot be 
boosted as much as in the case of a scale invariant spectrum, and the allowed range of parameters for our model (see 
Fig. [l^) is narrowed. 
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FIG. 9. Gravitational waves spectrum in a quintessence model with (IC2) initial conditions with a SUGRA-like potential. 
The spectrum has been normalized to be compatible with COBE at large scales. The spike in the evolution of the equation of 
state of the scalar field (see figure ^ yields to a loss [right] of about 20% in the amplitude of the spectrum at high frequency 
(i.e. > 10* Hz). 
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FIG. 10. Variation of the gravitational waves spectrum with the cut-off and with the epoch of transition between scalar field 
in kinetic regime and radiation [left] and contour plot [right] of the safe zone for nucleosynthesis for the parameter set (Tr, z,) 
(above the dotted and dot-dashed lines). Both plots are for a spectral index nr = and for maximum At allowed by COBE 
measurements. The solid line is obtained by imposing that the scalar field dominates at the end of infiation: points lying on or 
above this line have a corresponding solid power spectrum on the left plot. The horizontal dot-dash line is obtained by imposing 
that we are radiation dominated at nucleosynthesis (i.e. z* > 10^"). The diagonal dotted line is obtained by imposing that the 
energy density at the end of infiation is smaller than the Planck density. The diagonal dot-dashed line ( "GW background. . . " ) 
is obtained by imposing that f^ow < 10~^. 



VIII. CONCLUSION 



In this article, we have studied some properties of quintessence models with a non-minimally coupled scalar field 
among which the spectrum of gravitational waves. 

We have shown that such a quintessence field can behave as a fluid with ut < —1 and our models lead to —3 w < 
when the field dominates. We related the energy scale A of the potential to its slope a and to the scalar field energy 
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density today fi^. In particular, we showed that A is almost independent of 0° . The coincidence problem, i.e. the 
fact that fl^ ~ 1 implies a tuning of A (roughly the precision on A has to one order of magnitude higher than the 
one on $7^) which is however less severe than the fine tuning needed for a cosmological constant. This being fixed, 
the tracking mechanism allows to span a very wide range of initial conditions for the scalar field and there is no fine 
tuning in that respect. 

We then showed that the combined study of the gravitational waves energy spectrum and of their imprint on the 
CMB radiation temperature and polarisation enables to extract many complementary informations on the models: 

• the CMB mainly gives results on ^, and ut, 

• the energy spectrum gives results on the initial conditions of the scalar field. 

As pointed out in there is an excess of gravitational waves today if inflation ends by a kinetic phase. In that 

case, one has to check that both and fJcw are negligible at the time of nucleosynthesis and we relate the amount 
of gravitational waves today to the reheating temperature and the time of equality between the kinetic scalar era and 
the radiation era. 

We also pointed out that gravitational waves are damped by the anisotropic stress of radiation, which implies that 
the CMB anisotropy and polarisation spectra are lowered roughly by 10% for high multipoles. It was also shown 
that the amplitude of the gravitational waves spectrum for inverse power law potentials is 30% higher than for 
SUGRA-like potentials at high frequency. Indeed this is probably not detectable by coming experiments but it could 
ultimately lead to a signature of supergravity. 



ACKNOWLED GMENTS 



It is a pleasure to thank Pierre Binetruy, Nathalie Deruelle, Thibault Damour, David Langlois, Patrick Peter and 
Filipo Vernizzi for fruitful discussions. 



A.G. Riess et ai, Astronomical J. 116, 1009 (1998). 
P.M. Garvanich et ai, Astrophys. J. 509, 74 (1998). 
S. Perlmutter et ai, Nature (London) 3 91, 51 (1998). 



A. Melchiorri et ai, [ ELstro-ph/991144. 



Y. Mellier, Ann. Rev. Astron. Astrophys. 37, 127 (1999). 
R. Juszkiewicz et ai, Science 287, 109 (2000). 

N. Bahcall, J. P. Ostriker, S. Perlmutter and P.J. Steinhardt, Science 284, 1481 (1999). 

G. Efstathiou, Month. Not. R. Astron. Soc. 310, 842 (2000). 

S. Perlmutter, M.S. Turner and M. White, Phys. Rev. Lett. 83, 670 (1999). 



R.R. Caldwell, [astro-ph/9908168 



S. Weinberg, Gravitation and Cosmology (Addison Wiley and sons. New York, 1972). 
K. Coble, S. Dodelson and J.A. Frieman, Phys. Rev. D 55, 1851 (1997). 
R.R. CaldweU, R. Dave and P.J. Steinhardt, Phys. Rev. Lett. 80, 1582 (1998). 
P. Binetruy, Phys. Rev. D 60, 063502 (1999). 
P. Brax and J. Martin, Phys. Lett. B468, 4 (1999). 
P. Brax and J. Martin, [ ^stro-ph/991204e |. 



C. Wetterich, Astron. Astrophys. 301, 321 (1988). 
B. Ratra and P.J.E. Peebles, Phys. Rev. D 37, 3406 (1988). 
P.J. Steinhardt, L. Wang and I. Zlatev, Phys. Rev. D 59, 123504 (1999). 

G. Huey, L. Wang, R. Dave, R.R. Caldwell and P.J. Steinhardt, Phys. Rev D59, 063005 (1999). 
I. Zlatev, L. Wang and P.J. Steinhardt, Phys. Rev. Lett. 82, 896 (1999). 
J.-P. Uzan, Phys. Rev. D 59, 123510 (1999). 

N.D. Birrel and P.C.W. Davies, Quantum fields in curved space (Cambridge University Press, Cambridge, England, 1982). 
L.H. Ford, Phys. Rev. D 35, 2955 (1987). 
K. Maeda, Class. Quant. Grav 3, 233 (1986). 

F.S. Accetta, D.J. ZoUer and M.S. Turner, Phys. Rev. D 31, 3046 (1985). 
F. Lucchin, S. Matarrese and M.D. Pollock, Phys. Lett. B167, 163 (1986). 



20 



L. Amendola, Phys. Rev. D 60, 043501 (1999). 

S. Carol!, Phys. Rev. Lett. 81, 3067 (1998). 

T. Chiba, Phys. Rev. D 60, 083508 (1999). 

N. Bartolo and M. Pietroni, Phys. Rev. D 61, 023518 (2000). 

O. Bertolami and P.J. Martins, Phys. Rev. D 61, 064007 (2000). 

L.P. Grishchuk, JETP 40, 409 (1974). 

L.P. Grishchuk, Sov. Phys. Usp. 31, 940 (1988). 

B. Allen, in Relativistic gravitation and gravititional radiation, edited by J. -A. Marck and J. -P. Lasota (Cambridge Uni- 
versity Press, Cambridge, England, 1997). 

A. Vilenkin and E.P.S. Shellard, Cosmic strings and other topological defects (Cambridge University Press, Cambridge, 
England, 1994). 

A. Abramovici et ai, Science 256, 325 (1992). 

C. Bradaschia et ai, Nucl. Instrum. Meth. A289, 518 (1990). 

Y.R. Jafry, J. Cornelisse and R. Reinhard, ESA Journal 18, 219 (1994). 
V. Rubakov, M. Sazhin and A. Veryaskin, Phys. Lett. B115, 189 (1982). 

A. A. Starobinsky, Sov. Astron. Lett. 9, 302 (1983). 

L.F. Abbott and M.B. Wise, Astrophys. J. 282, L47 (1984). 
M. White, Phys. Rev. D 46, 4198 (1992). 

M.S. Turner, M. White and J.E. Lidsey, Phys. Rev. D 48, 4613 (1993). 

W. Hu and M. White, New Astron. 2, 323 (1997). 

M. Giovannini, Phys. Rev. D 60, 123511 (1999). 

M. Giovannini, Class. Quant. Grav. 16, 2905 (1999). 

P.J.E. Peebles and A. Vilenkin, Phys. Rev. D 59, 063505 (1999). 

B. Spokoiny, Phys. Lett. B315, 40 (1993). 

M. Giovannini, Phys. Rev. D 59, 121301 (1999). 

A. Me lchiorri, M.V. Sazhin, V.V. Sh ulga and N. Vittorio, Astrophys. J. 518, 562 (1999). 
MAP, tittp : //map .Rsf c .nasa.gov/ 



Planck, [http : //astro . estec . esa.nl/Planck/ 



H. Kodama and M. Sasaki, Prog. Theor. Phys. Suppl. 78, 1 (1984). 
J.-P. Uzan, Class. Quant. Grav. 15, 1063 (1998). 
L.H. Ford and L. Parker, Phys. Rev. D 16, 1601 (1977). 

C.W. Misner, K.S. Thorne and J. A. Wheeler, Gravitation (Freeman, San Francisco, 1973). 

W. Hu and M. White, Phys. Rev. D 56, 596 (1997). 

R.K. Sachs and A.M. Wolfe, Astrophys. J. 147, 73 (1967). 

E. Newman and R. Penrose, J. Math. Phys. 7, 863 (1966). 

J.N. Goldberg et al, J. Math. Phys. 8, 2155 (1967). 

K.S. Thorne, Rev. Mod. Phys. 52, 299 (1980). 

P.J.E. Peebles, Astrophys. J. 153, 1 (1968). 

S. Seager, D.D. Sasselov and D. Scott, Astrophys. J. 523, LI (1999). 
C.-P. Ma and E. Bertschinger, Astrophys. J. 455, 7 (1995). 
G.T. Gillies, Rep. Prog. Phys. 60, 151 (1997). 
T. Damour and F. Dyson, Nucl. Phys. B480, 37 (1996). 
T. Damour and G. Esposito-Farese, Class. Quant. Grav. 9, 2093 (1992). 
P.G. Ferreira and M. Joyce, Phys. Rev. D 58, 023503 (1997). 



M. Tegmark ar id M. Zaldarriaga, Astrophys. J. , to appear (2000), [ ^stro-ph/0002091 



H. Ziaeepour, [ |astro-ph/0002400[ . 
T. Chiba, T. Okabe and M. Yamaguchi, [ ELstro-ph/9912463| . 



A. Dimitropoulos, Astron. Astrophys. Trans. 18, 543 (1999), | |gr-qc/9805016 
J. Silk, Astrophys. J 151, 459 (1968). 

W. Hu, PhD thesis. University of California at Berkeley (1995). 
L.P. Grishchuk, Phys. Rev. D 48, 3513 (1993). 
M. Joyce, Phys. Rev. D 55, 1875 (1997) 
L.P. Grishchuk, | lgr-qc/0002035 



21 



